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Abstract. We give bounds on the rate of convergence to equilibrium of the symmetric 
simple exclusion process in Z d . Our results include the existent results in the literature. 
We get better bounds and larger class of initial states via a unified approach. The method 
includes a comparison of the evolution of n interacting particles with n independent ones 
along the whole time trajectory. 



1. Introduction. 

The rate of convergence to equilibrium is one of the main problems in the theory of 
Markov processes. It has recently attracted the attention of many authors in the context 
of symmetric conservative particle systems in finite and infinite volume. In finite volume 
the techniques used to obtain the rate of convergence to equilibrium rely mostly on the 
estimation of the spectral gap of the generator. In general, one shows that the generator 
of the particle system restricted to a cube of length iV has a gap of order N~ 2 in any 
dimension. This estimate together with standard spectral arguments permits to prove 
that the particle system restricted to a cube of size iV decays to equilibrium in L 2 at 
the exponential rate exp{— ct/N 2 }. This approach has been successfully extended to 
the infinite volume setting and permitted to prove L 2 -polynomial decay to equilibrium. 
The method, however, does not give any information on the rate at which the system 
converges to equilibrium starting from an arbitrary configuration or from arbitrary initial 
distributions. In this article we consider the symmetric simple exclusion process in which 
some explicit computations can be performed. 

The symmetric simple exclusion process was introduced by Spitzer (1970). Informally 
one can describe the process following the so called stirring representation. Fix a sym- 
metric transition probability p on Z d : p(x,y) = p(y,x), p(x,y) > 0, ^2 v p(x,y) = 1 
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for all x in Z d . We assume that the transition probability is translation invariant, 
p{x,y) = p(0,y — x) and that it is indecomposable in the sense that for each z in Z d , 
there exists an integer n and a sequence = zq, . . . , z n = z such that p(zi, Zi+i) > for 
alH = 0, . . . , n — 1. For each pair of sites x, y such that p(x, y) > 0, consider a Poisson 
point process, denoted by N Xjy (t), with rate p(x, y). Assume that these processes are all 
independent. As initial state, fix some configuration rj with at most one particle per site. 
Thus, the configuration rj is a collection rj = {rj(x), x G Z}, where rj(x) = 1 indicates that 
the site x is occupied for the configuration rj and rj(x) = means that the site is empty. 
To obtain the state of the process from the initial states and from the Poisson point 
processes, we proceed as follows. Each time the Poisson process N XiV increases by 1, we 
interchange the variables rj(x) and rj(y). Notice that if both site x and y are occupied or 
if both are vacant before the jump of N XjV , the configuration remains unchanged after 
the jump. In the other two cases the modification can be interpreted as the jump of 
a particle from the occupied site to the unoccupied site. Notice also that at each later 
time each site is occupied by at most one particle and that the total number of particles 
is conserved by the dynamics. This explains why this stochastic dynamics is called the 
symmetric simple exclusion process and in which sense it is conservative. 

For p E [0, 1], denote by v p the Bernoulli product measure with density p. This is the 
probability measure on the configuration space {0, 1} Z obtained by putting a particle at 
each site with probability p independently of the other sites. Liggett proved that all the 
invariant measures for the symmetric simple exclusion process are convex combinations 
of the Bernoulli measures v p (cf. Liggett (1985)). 

When one tries to study convergence to equilibrium for arbitrary initial configurations, 
the first issue is the very question: "what to prove?" The first attempt is to try to prove 
the existence of a function h(t) that decreases to as t | oo and the existence of a norm 
I /I such that for all initial configuration rj chosen according to the invariant measure v p 
and for all cylinder function / 

IW-"p/I < Ill/I *(*)■ 

Here {S(t),t > 0} stands for the semigroup of the symmetric simple exclusion process 
and Srj for the probability measure concentrated on rj so that 5 v S(t) is the distribution of 
the process at time t starting from rj. Also, for a bounded function / and a probability 
measure p, pf stands for the expectation of / with respect to p. 

This is clearly not true as for any fixed e > and any fixed t, one can always choose 
a set of configurations with v p positive probability for which 

\8 v S(t)[rj(0)]-p\ > e. 

Indeed, it is enough to consider configurations rj whose sites in a cube of length t around 
the origin are all occupied. In general there is no hope to have uniform almost sure 
convergence in conservative systems. In spin flip systems the equilibrium is attained 
locally in an independent way for distant regions. For this reason one can hope to get 
almost sure convergence to equilibrium in spin-flip systems. 
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The second attempt is to fix a configuration rj, to choose the density p depending on 
the initial configuration: pi = SrjS(t)[ri(0)] and to compute 

\5 n S{t)f-v pl f\. (1.1) 

This approach cannot give a bound better than c(/)t -1 / 2 in any dimension. This is not 
satisfactory because, in view of the decay to equilibrium in L 2 , one expects to obtain 
estimates of order t~ d l 2 in dimension d. To check that this formulation can not give 
bounds better than t -1 / 2 , consider the configuration r\ such that i](x) = 1 if and only 
if x\ > 1. Here x\ stands for the first coordinate of x. In this standard duality 

argument (that will be explained in section 2) gives that pi = P[X^ > 1], where X° is 
a one dimensional symmetric random walk that starts from the origin. Fix the cylinder 
function f(rj) = ??(ei), where e^, 1 < i < d, is the canonical basis of R d . By duality and 
translation invariance, 5 v S(t)f = P[X® > 0]. In particular, for this cylinder function, 
(1.1) is equal to P[X® = 0], which is of order t -1 / 2 . 

At this point we have two possibilities. We may of course impose some regularity 
conditions on the initial configuration r\ (assume for instance that it is periodic) or to 
average the difference \S v S(t)f — u p f\ (or a power of this difference) with respect to some 
measure v that has asymptotic density p and some nice correlation properties. The other 
possibility is to take advantage that we are in the context of exclusion process, where all 
cylinder functions can be written as linear combinations of functions of type n^e/i v( x ) 
for some finite set A. In this CclSG, cl natural quantity to investigate is 

xEA 

where pl(x) = 5 v S(t)r](x). These are the so-called ^-functions introduced by Ferrari, 
Presutti, Scacciatelli and Vares (1991). A related quantity, also very natural in the 
context of the symmetric simple exclusion process is the difference 

5 v S(t) J] V (x) - n S v S(t) V (x) . (1.2) 

xeA xeA 

We shall follow the two directions just mentioned. We shall first prove that the 
difference (1.2) can be expressed in terms of quantities related to independent random 
walks with transition rate p(x, y). This is the content of Theorem 2.2. This fact together 
with some elementary bounds on the transition probability of symmetric random walks 
will permit to obtain sharp estimates on the integral of powers of (1.1) with respect to 
translation-invariant measures that have density p and polynomial decaying correlations. 
This is the content of Theorem 2.1. 

The comparison between n random walks interacting by exclusion and n independent 
random walks was studied by Bertein and Galves (1977), De Masi, Ianiro and Presutti 
(1982), De Masi and Presutti (1983), De Masi, Ianiro, Pellegrinotti and Presutti (1984), 
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Ferrari and Goldstein (1988), Ferrari, Presutti, Scacciatelli and Vares (1991) and Andjel 
(1994). We give a short and unified presentation which includes most of the above results. 

A simple exclusion process in a finite box can be understood as a random walk in 
a finite set. Convergence to equilibrium for finite-state Markov processes have received 
new attention lately. We quote for instance Aldous (1983), Diaconis (1988), Diaconis 
and Stroock (1991). For the symmetric simple exclusion in a finite box, Quastel (1992) 
computed the spectral gap of the generator, which gives the L 2 exponential rate of 
convergence to equilibrium in finite volume. Using this result, Bertini and Zegarlinski 
(1998) proved polynomial L 2 -convergence to equilibrium (time-correlation decay for the 
system in equilibrium) in infinite volume. Janvresse, Landim, Quastel and Yau (1999) 
proved a analogous result for the symmetric zero range process. 

Cancrini and Galves (1995) obtained an upper bound for the rate of convergence of 
symmetric simple exclusion processes starting either from a periodic configuration or 
from a stationary measure satisfying mixing conditions. This result was extended to 
the one-dimensional nearest-neighbor zero-range process with rate g(k) = l{k > 1} by 
Galves and Guiol (1997). 

2. Notation and results. 

Let p(x, y) be a symmetric, translation-invariant, irreducible, non-negative real matrix 
on Z d such that 

y 

To each pair of sites (x, y) of Z d attach a Poisson process of rate p(x, y) denoted by 
N x ,y{t). We adopt the following convention: dN x ^ y (t) = 1 when there is an event of the 
corresponding Poisson process at time t. As a function of these Poisson processes, for 
s < t, we define a family of bijections of Z d in itself, £|: Z d — > Z d , such that 

C(x)=x and d^ = J2mr y -^ t }dN s x Jt), 

x,y 

where for any bijection £: Z d — > Z d , 

{i(z) if z^x,y, 
£(v) Xz = x, 
£{x) if z = y 

and TV* (£) = N X}V (s + 1) — N XjV (s). £ s is thus the identity at time s and if there is a 
mark of the Poisson process N XyV at time s + t, £ s (x) and $, s (y) change their value in the 
following way : (x) = £%_(y), Ctiu) = £t-( x )- ^ n particular, for any fixed x, the process 
£, s (x) is a stochastic process that may have very long jumps even if the rates •) are 
near est- neighb or . 

The simple exclusion process with initial configuration rj is defined as 
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for all t > 0. Denote S(t)f(rj) = E/(7#). Here and below P and E stand for the 
probability and the expectation with respect to the Poisson point processes. The set of 
extremal invariant measures for this process is the set of Bernoulli product measures v p 
indexed by the density p G [0, 1] . The main results of this article concern the convergence 
to those measures from different initial conditions. They are stated in the next theorem. 

Fix a > 0. We say that a measure v on {0, 1} Z has a-decaying correlations if there 
exist a finite constant C such that 

E v [fg]-E v [f]E v [g) < C||/ Ukj" ' 



d(A /} A ff )« 

for every cylinder functions /, g. In this formula, Af,A g stands for the support of the 
cylinder functions /, g and d(A, B) for the distance between two subsets A, B of Z d . 

Theorem 2.1 . Let v be a translation-invariant probability measure in {0, 1} Z with 
density p (that is E v [r)(x)] = p) and with a-decaying correlations for some a > d. Then, 
for each cylinder function f there exists a constant c(f) such that for all t > 0: 

(a) Weak convergence: 

\uS(t)f-u p f\ < c(f)(l+t)- d / 2 

(b) L p convergence for 1 < p < 2: 

J v(d V )\8 v S(t)f-v p f\v < C (f)(l + t)- d ^ 4 . 



(c) Lp convergence for p > 2: 

J v(d V ) \S v S(t)f - u p f\* < c(f) (1 + t)~ d ' 2 . 

Remarks. Part (a) in the theorem above improves the bound obtained by Cancrini 
and Galves (1995) in two directions: on the one hand, we only ask for polynomial decay 
of correlations in the initial measure, while [CG] requires exponential decay. On the 
other hand, in [CG] the upper bound includes a (spurious) logt in the numerator. Part 
(b) generalizes a result obtained by Bertini and Zegarlinski (1996) in two senses: [BZ] 
proved (b) for p = 2 with product initial measure, while part (b) requires only the initial 
measure to have polynomial decay of correlations, which is of course satisfied by v p , and 
part (b) is proved for 1 < p < 2. Finally, since v p is an equilibrium measure, part (b) 
can be read as the time-decay of correlations of the system in equilibrium. 

The proof of Theorem 2.1 is based on the self-duality of the symmetric simple exclusion 
process that we now explain. For < s < t, let X%'* denote the position at time t — s 
of the particle sitting at x at time t : X®'* = Cl~ s ( x )- m contrast with ^ s (x), X x,t is a 
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random walk with transition rate p(-, •). Moreover, for any subset A of Z d , {X x,t , x G A} 
evolves as symmetric exclusion random walks and 

P[t/?(x) = 1, x G A] = Pfaprf'*) = 1, x G A] 

for every finite subset A. This is the so-called self-duality relation of the symmetric 
simple exclusion process. 

In the following theorem we compare the evolution of n particles interacting by ex- 
clusion with n independent particles. The proof is inspired in a similar result by Ferrari 
and Goldstein (1988), where exclusion processes with birth and deaths were considered, 
and a result by Ferrari, Presutti, Scacciatelli and Vares (1991). This can be seen as 
a probabilistic version of the "integration by parts formula". See Proposition 8.1.7 of 
Liggett (1985) and display (1) of Andjel (1994), for instance. 

Theorem 2.2 . For any vector x = (x\, . . . , x n ), 

n n 

En^)-n E ^) 

i=i i=i 

for all t and any configuration r\. In the above formula the summation is carried over all 
V = (j/i, • • • , Vn) such that y { ^ y d for i ^ j and 

p?(y) = E^(y). 

In the next theorem we apply the previous result. 

Theorem 2.3 . Assume p(-, ■) has a finite second moment. 

(a) There exists a finite constant C such that for all t > 0, 

n n 

su P |e n^(^) - n Er ^)| ^ cr ^) > 
v i=i i=i 

where R^t) = log(l + t)/y/T+t, R 2 (t) = log(l + t)/(l + t) and R d (t) = (1 +t)~ 1 
for d>3. 

(b) For each positive function <p : R + — > M, + let 



{77: forallx,yeZ d , \x - y\- 2 (p v t (x) - p v t (y)) 2 < <p(t)} 
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Then, for each decreasing function ip such that ip(t) < Cq(1 + t)~( li + 2 )/ 2 for some 
finite constant Co, there exists a finite constant C\ such that for all t > 



sup ll Y[v?(xi)-]1®V?( 



.1, i 



i=l 



i=l 



<Ci 



(i + ty/ 2 ' 



(c) Let v satisfy the conditions of Theorem 2.1. Then there exists a finite constant C 
such that for all t > 



i=l 



i=l 



< c 



(l + t) d / 2 



Remarks. Part (a) was obtained by Ferrari, Presutti, Scacciatelli and Vares (1988). 
Let r\ be a periodic configuration on Z d . It is easy to show that r\ belongs to X<^ for some 
ip satisfying the assumptions of Theorem 2.3 (b). Part (b) includes a result in Landim 
(1999), where the case of finite initial r\ is considered. 

3. Proofs. 

Before proving the theorems, we state an estimate that will be needed several times 
in the sequel. It is based on the classical negative correlations property of the symmetric 
exclusion process. 

Lemma 3.1 . Let $ : Z d \ {0} — > M + be a summable function. Then, there exists a 
constant C such that for all set of different sites {x\, . . . , x n } and for all t > 0, 

E E v[{x? 1 '\---,x? n ' t } = {vu---,v n }]*(y j -Vk) < C{i + ty d '\ (3.1) 

jj^k yi,...,y n 

where the second sum is carried over the set {j/i, . . . , y n : y a 7^ yt>, for a ^b}. 

Proof. If we fix j ^ and sum over yi for % 7^ j , k, we obtain that the left hand side of 
(3.1) is equal to 

E E n{ X t ut ,---,X^ t }D{y j ,y k }}^( yj - yk ). (3.2) 
The previous probability can be decomposed as 

E p [W B,t ^t X6 ''} = {yi,y*}]- 
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By the correlation inequality between symmetric exclusion random walks and symmetric 
independent random walks (cf. [L], Theorem VIII. 1.7), each of theses probabilities is 
bounded above by 

P[Xf-'*G{y i ,y fc }]P[Xf"'G{y i ,y*}] 

= {Pt(Vj -x a ) +Pt(Vk -x a )}{p t (yj -x b ) +pt(yk -Xb)} 

provided pt(x, y) stands for the probability of a continuous-time random walk with tran- 
sition probability p(-, ■) to be at y at time t if it starts from x at time and provided 
Pt(x) stands for p t (0,x). The sum (3.2) is thus bounded above by 

EE E ®(yj ~Vk) {Pt(yj ~ x a ) +Pt(yk - x a )}{p t (yj - x b ) +p t (y k - x b )} . 

Since $ is summable and since pt(x) is bounded above by C/(l + t) d ^ 2 for some finite 
constant that depends only onp(-), this expression is bounded above by C(n, p)(l+t)~ d ' 2 , 
which concludes the proof of the lemma. □ 

Proof of Theorem 2.1 (a) Let v be a translation-invariant probability measure that 
has a-decaying correlations. By duality, for any n > 1 and any distinct sites x±, . . . , x n , 

^u[vt(xi) = 1, % = 1, ...,n] - p r 

= | E V[{Xr'\...,X^ t } = {y 1 ,...,y n }}{v{ V (y l ) = l^ = l,...,n}-p n }\. 
yi,---,y n 

Since v is translation-invariant with mean p and has a-decaying correlations, there exists 
a finite constant C such that 

V{rj( yi ) = l,i=l,...,n}-P n \ < Cj2 r y * Ja • 
The right hand side of the previous expression is thus bounded above by 

C H E n{xr>\...,x^ t } = {yu...,y n }]T—^- 

j^k 91--* |yj yfc| 

To conclude the proof it remains to apply Lemma 3.1 for the function $(y) = |y| _a - We 
can do this because in the above expression yj ^ yk- □ 

Proof of Theorem 2.2. Fix n distinct sites x\, . . . , x n and define {Y^J' S : % = 1, . . . , n} 
as a family of independent random processes with the same marginal distribution as 
Xf i,s respectively. 



CONVERGENCE TO EQUILIBRIUM OF THE EXCLUSION PROCESS 



9 



We realize the motion of the Y process by considering a family of Poisson marks as 
the one defined by N XiV (t), called N® y (t), independent of the precedent one as follows. 
Processes X uses only the N marks. Process Y uses the N marks if only one Y particle is 
concerned by the jump. This means that if there is a Y particle at x and no Y particle at 
y at time t, then the process uses the marks of N Xjy (t) and ignores the marks of N® y (t). 
If the jump concerns two Y particles, say particles % and j with % < j, then the Yi particle 
uses the N marks and the Yj particle uses the N° marks to jump over the position of the 
Yi particle and the N marks to jump to any other position. Thus if particle Yi is at x 
and particle Yj at time t with % < j, particle Yi uses the marks of N Xiy (t) while particle 
Yj uses the marks of (t). 

We just gave a coupling between a system of n exclusion and n independent particles. 
This means that we realized the two processes in the same probability space (the one 
generated by the product of the Poisson processes {N Xyy {t), N® y (t), ; x,y e Z d , t > 0}) 
in such a way that the marginal distributions are those desired for both processes. We 
continue using P and E for the probability and expectation with respect to the product 
of the Poisson processes. 

By definition of the symmetric simple exclusion process, 

n n n n 



Let Ti be the first instant that two Y particles occupy the same site. We say that a 
collision occurred at that time. Before the collision each X particle occupies the same 
place of the corresponding Y particle. For s < t let 

n n 

i(x u . . . , x n , s, t) = n nix?*) - n ntrg'*) ■ 

i=l i=l 

We want to compute the expectation of I(xi, . . . , x n , t) = I(x\, . . . , x n , t,t). If t < Ti, 
I(xi, . . . , x n , t) is zero because the trajectories of the X and Y process coincide. On the 

other hand, since for all < s < t, Y?j '* = Y^ s ' t_S , X? '* = l£''''" s , on the set 
{Ti<t}, 

n x ■ t n x ■ t 

i( X1 ,...,z n ,t) = Un(xli ,t " Ti )-n^54 , '" Tl ) • ( 3 - 3 ) 
i=i i=i 

( 1 if the collision at T\ occurs due to a N mark 
[ if the collision at T\ occurs due to a iV mark 

1 < i < j < n} stand for the labels of the particles involved in the 



Let 

Ex = 

and let Z\ G 
collision at time Ti : 



Zi = if and only if = Y*£ and i < j 
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Assume that the first collision is due to a collision between particles % and j with i < j. 
In this case, at time Ti, the position of particles Y£ k,t and X Xk ' t coincide for k ^ i, j. 
Moreover, if the collision occurred due to a N mark, i.e., due to a jump of particle Yi 
over particle Y h Y?l' 1 = X- 3 T = Y" i . In the case where the collision occurred due to a 

J * ; -*- 1 j ; 1 Ji-'l 

N° mark, i.e., due to a jump of the Yj particle over Y{, a similar identity holds with the 
roles of i and j exchanged. In particular, on the set {Ti < £}, 



_ -y x k>t x rp -V^^fe;* . rp 

rNnA; " x ) = = n ^a; " x ) >< 

x{^=owAr m^a; ) 

X ■ , t X A , t 

XrJ ,t-Tl XrJ ,t-Tl "I 

+ i{^ = i}^A m^a )}• 



^X T ^ ,t—Ti, 



We now add and subtract in the right hand side rj(Y i t _ 1 Ti )v(Yj t-T t ) to recover 



rii<A:<n ^(^fc t-Ti )• After this step we obtain that on the set {Ti < £}, 



n ^ A ) 

fc=l 



k=i 



i<j 



x (l^ = 0} [tfy, A )r / (y J . A ) - ^ A )^ y ,- A ) 



X£'*,t-Ti 



X^^t-T, 



Xrj? , t Tl . . Xrj? . t T\ XrJ ' , t T\ X rJ , t 7\ 

+ 1{E, = 1} „(K W A ) - »y(y 4 A, W. 



x i — t 1 ' — ^ 

Denote the right hand side of this expression by g(Ei, Ti, Z\,Xt£ , Y_ t _j, ). It follows 
from this identity and (3.3) that 

I(x,t) = l{Ti KtyiUL^fit-T!) - l{T 1 <t}g(E 1 ,T 1 ,Z 1 ,X%,Y^*~ Tl ). 



For a positive integer £, define Ge = cr{T m , Z m , X^y 1 , . . . , J^'* : m = 1, . . . , £}. Since 
Ti is a Bernoulli random variable with parameter 1/2 independent of C?i and since the 
y particles evolve independently, by the strong Markov property on the set l{Ti < t} 
we get that 

E(g(E 1 ,T 1 ,Z 1 ,X^ 1 Yf^; t ~ Tl )\g 1 ) = g{T u Z lt . . . , X*-'*) , 
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where 

g(u,^j},x 1 ,...,x n )= 1 - [E HY^ u ))-E( V (Y^ u ))] 2 J] E (r^j) . 

From the two previous identities, we obtain that 

E(I(x lt ... t x nt t)\g 1 ) = l{T 1 <t}I(X%' t ,... t X%> t ,t-T 1 ) 

- l{Ti < Zx, . . . , X*f) . 

Repeating the argument for I > 2, we get the following expression for the expectation 
of/: 

M(t) 

E (J(xi, ...,i m t)) = -E^ </(T £ , Z,, . . . , X*;'*) , 

where T^E^Z^ are defined as Ti, Ei, Zi inductively and M(£) = X^>i — ^ s 
the number of collisions occurred by time t. Notice that since g is positive, the above 
expression implies immediately that the distribution of rj t is a measure with negative 
correlations for any initial rj. 

For i < j, denote by T\' 3 the instant of the £-th collision of particles i,j and by M l > J (t) 
the number of collisions up to time t of particles i,j so that M{t) = J2i<j M l ' J (t) and 

M*' j (t) 

E (I( X1 , . . . , x n , t)) = - E E E ^ - V r''- • • • » » ( 3 - 4 ) 



where 



9 hJ (^y) = \ [e M^i'J) -e(^j)] 2 n E ^i) • ( 3 - 5 ) 



For i < j, M l ' J (t) is a Poisson process with rate 2 = y}p(Vi, Dj)- Moreover, 

for any function h(u,x), 

E h(r e > 3 ,x%t) = / ^.a'l^w- 

In particular, taking expectations, 

M^(t) t 

E E //'-'(Vr-'.^'J = 2/ ^P(w^i) p (^f = v)^ j '(a»y)*'- 
This together with (3.4), (3.5) concludes the proof of the theorem. □ 
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Proof of Theorem 2.3. To show item (a) we use Theorem 2.2. Since we are taking 
supremum over rj, we can cancel the product Ylk^i j Pt-siVk)- We get the following upper 
bound 

Ef E p(y i ,y j )F({X^,...,X^} = {y 1 ,...,y n })x 

1<J J0 yi ,...,y n (3.6) 

x {pl-siVi) - Pl-siVi)) 2 ds 

Denote by p t (x,y) the probability for a random walk with transition probability p(-, •) 
starting at x to be at y a time t. By Lemma 3.2 below, there exists a universal constant 
C\ such that 



Ip2(2/«) — P2(2/j)l = I 5E ^ (Ps(Vii z ) -Pa(Vj,z)) 

zei d 

for all s > 0. By Lemma 3.1 with $(y) = p(0,y)\y\ 2 , which is summable by hypothesis, 
(3.6) is bounded above by 

C [\l + s )- d / 2 — -± ds < CR d (t). 

Jo 1 + t-s 

Here and below C stands for a finite constant that may change from line to line. This 
shows item (a). To show (b), again we can cancel the product Ylk^i j Pt-siVk) an d need 
to bound (3.6) for rj in X^,. For those 77's we have 

(Pt-siVi) -Pt-siVj)) 2 < <P(t- s)\yi-yj\ 2 . 
Applying Lemma 3.1 with <E>(y) = p(0, y)\y\ 2 , we bound (3.6) by 

C [ (1 + s)~ d/2 (p(t - s) ds < C(l + t)~ d/2 
Jo 

because (p(t) was assumed to be bounded above by C(l + t)~( d + 2 )/ 2 . This proves (b). 

To show (c), as before, we need to compute the expectation with respect to v of 
(3.6). Hence using Fubini we need to compute v{pl_ s {yi) — Pt-siUj)) 2 - Since pi (y) = 
Ylz v( z )Pt(y, z ), developing the square we get 

v(pt- s (yi) - Ptsivj)) 2 

= u (v( z i)v(z 2 )) (pt(yi,zi) -Pt(yj,zi))(p t (yi,z 2 ) -p t (yj,z 2 )) 

Zl,Z2 

= P^2(pt(yz,z) -pt{y^z)f 

z 

+ $E ( u v(zi)v( z 2) - P 2 ) (pt(yi, zi) - Pt(Vj, zi)) (pt(yi, z 2 ) - p t (yj,z 2 )) 

Z1t^Z 2 

+ 5E p 2 (p*(f<> Zi ) ~ pt(yji z i)) (pt(yi, z 2 ) - pt(yj, z 2 )) . 

zi=£z 2 
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It follows form Lemma 3.2 below that the sum of the first and third line of the previous 
expression is bounded above by Cp(l — p)\yi — yj\ 2 (l + t)~^ d+1 - ) / 2 , while the second, by 
Schwarz inequality, is bounded above by 

E \ u v( z i)v(z2) - p 2 \{pt(yt,z 1 ) -pt(yj,z!)} 2 . 

Since v has a-decaying correlations, by Lemma 3.2, this expression is bounded above by 
C\yi - yj\ 2 {l +t)" (d+1)/2 . Therefore, (3.6) is less than or equal to 

c (p)Y,f E p^y^-yj\ 2 nKf t = yj (1+t _ 1 s)(d+1)/2 d*. 

i<3 yi,---,Vn ^ ' 

It remains to apply Lemma 3.1 to the function $(y) = p(0 7 y)\y\ 2 to obtain that the 
previous expression is bounded above by 

rt 1 1 



Jo + s 



ds . 



s )d/2 ( 1+t _ s )(d+l)/2 

This concludes the proof of the theorem in dimension d > 2. In dimension 1, we need 
(3.7) instead of Lemma 3.2 to prove the estimate. □ 

We turn now to Theorem 2.1. 

Proof of Theorem 2.1. We already proved part (a). In order to prove (b), fix 1 < p < 2. 
It is enough to consider cylinder functions of type /a = rix6A^( ;r ) f° r a fi n he subset A 
of Z d . For such a cylinder function, 5 v S(t)fA = E ELe^i Vt( x ) an d ^ P fA = provided 
\A\ stands for the number of sites of A. Since < /a < 1 and since a p < a for < a < 1, 
by Theorem 2.3 (c), 



J Hdv)\S v S(t)f A -is p f A \P < 2P- 1 J u(d V ) \ Y[E V ?(x)-p\ A \\ P + C(l + t) 



-d/2 



x€A 

for some finite constant C. Introducing intermediary terms, the first term on the right 
hand side is bounded above by 

p 



C(A)J2 J v(dr,)\Erfl(x)-p 



xeA ' 

Since p < 2, By Holder inequality, the previous expression is bounded above by 

2,p/2 



xeA J 



To show that this expression is bounded above by C(A)(l+t)~ dp ^ 4 it is enough to expand 
the square and to recall that v has a-decaying correlations and that pt(x) is bounded 
above by Ct~ d l 2 uniformly in x. The proof of (c) is exactly the same. This concludes 
the proof of the theorem. □ 

We conclude this section with an estimate on the transition probability of symmetric 
random walks. 
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Lemma 3.2 . Let {X t ,t > 0} be a random walk on Z d with transition probability p(-, ■) 
satisfying the assumptions stated in the beginning of the article. Then, there exists a 
finite constant C such that 

J2 \pt(0,x + e z ) -p t (0,x)\ < C(l + t)~ 1/2 and 

|pt(0,x + ei )-pt(0,x)| 2 < C(l + t)- (1+d)/2 . 

xei, d 

Here {ei, 1 < i < d} stands for the canonical basis ofM. d . 

Proof. This result can be proved by a coupling argument indicated to us by E. Andjel 
or through the local central limit theorem. With slightly stronger assumptions on the 
moments of p(-), the local central limit theorem gives better estimates (of type t~( d + 2 )/ 2 ) 
for the second term. In dimension 1, with the assumptions of indecomposability and finite 
second moments, we have 

J2\Pt(0,x + l)-p t (0,x)\ 2 < C(l+t)-^ (3.7) 

for some 5 > 0. 

The coupling argument is as follows. Consider two random walks X}, X 2 with tran- 
sition probability p(-) on Z d and with initial states and e\. Assume that p(ei) > 0. 
We couple these two random walks in the following way. If the process (Xl,Xf) is at 
(x 1 ,^ 2 ), for y 7^ ei it jumps to (x 1 + y, x 2 + y) at rate p(y), it jumps to (x 1 + ei,x 2 ) 
at rate p(ei) and it jumps to (x 1 ,^ 2 + ei) at rate p(ei). We proceed in this way until 
they meet. From this time on, they jump together. With this coupling the difference 
X\ — X 2 is a nearest-neighbor, symmetric, one-dimensional random walk with absorption 
at the origin that starts from —1. A well known bound gives that the probability for 
this one-dimensional random walk to have not reached the origin before time t decays 
as t~ x l 2 . This estimate permits to prove the lemma in the case p(ei) > 0. In the other 
case, since the transition probability is assumed to be indecomposable, there exists an 
integer n and a sequence = xq, x\, . . . , x n = e\ such that p{xi, Xi+\) > and we may 
proceed in a similar way. Notice that this proof does not require any assumption on the 
moments of p. □ 
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